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Abstract 



Canonization of F{R) theory of gravity to explore Noether symmetry is performed treating R^G{^ + ^ + -^) = 
■ as a constraint of the theory in Robertson- Walker space-time, which implies that R is taken as an auxiliary 

variable. Although it yields correct field equations, Noether symmetry does not allow linear term in the action, 
and as such does not produce a viable cosmological model. Here, we show that this technique of exploring 
Noether symmetry does not allow even a non-linear form of F{R) , if the configuration space is enlarged by 
O [ including a scalar field in addition, or taking anisotropic models into account. Surprisingly enough, it does 

not reproduce the symmetry that already exists in the literature [T5] for scalar tensor theory of gravity in the 
~-i ■ presence of R' term. Thus, R can not be treated as an auxiliary variable and hence Noether symmetry of 

arbitrary form of F{R) theory of gravity remains obscure. However, there exists in general, a conserved current 
for F{R) theory of gravity in the presence of a non-minimally coupled scalar-tensor theory [261I27| . Here, we 
briefly expatiate the non-Noether conserved current and cite an example to reveal its importance in finding 
^ ■ cosmological solution for such an action, taking F(R) oc R^ . 



PACS 04.50. -l-h 



0^. 

m; 

■ 1 Introduction 

o: . . . 

fSl . Increasing interest in F(R) theory of gravity (see [Tj for a recent review) initiates to probe deeply into it and 
' to explore its merits-demerits from different angles. Cosmological consequences of F{R) theory of gravity have 
been explored |2j taking into account different sets of field equations which appear through the metric variation 
and Palatini variation formalisms. Although, both set of field equations lead to late time cosmic acceleration 
without the requirement of dark energy in the form of scalar fields, different results have emerged following these 
' two different techniques. For example, unification of early inflation and late time acceleration [3] on one hand, 
and violent instability [4], cosmologically non- viability [5], big bang nucleosynthesis and fifth- force constraints 
altogether [6], on the other, are the outcome of metric formalism. On the contrary, no instability appears as 
such in Palatini formalism [7. and it leads to correct Newtonian limit [5], while curvature corrections are found to 
induce effective pressure gradients which create problem in the formation of large-scale structure [9]. Further, even 
density perturbation in the two techniques, produce different results jlOI . Nevertheless, it has been shown that 
under conformal transformation the two formalisms are the same dynamically [IV and also are equivalent for a 
large class of theories [12] . Despite such contrasting results on one hand and equivalence of the two formalisms on 
the other, at the end it is required to choose some particular form of F{R) out of indefinitely many. This is only 
possible by imposing Noether symmetry. In the present work we discuss this issue in the metric variation formalism. 



Noether symmetry when applied for the first time in scalar-tensor theory of gravity [13j . the idea was to find a 
form of the potential that might give rise to a cyclic co-ordinate and hence a conserved current. That is, out of 
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indefinitely many choice of tlie potential, Noether symmetry selects one, corresponding to which there exists a 
cyclic coordinate that can't be found a-priori through inspection. In the Robertson- Walker background metric 
the exponential form of the potential thus found was rather encouraging, because it could trigger inflation in 
the early Universe. Since then, it has been applied in nonminimally coupled scalar-tensor theory of gravity both 
in the background of isotropic and several anisotropic models (14j . together with nonminimally coupled scalar- 
tensor theory in the presence of term in the action PJ] and recently, to classify different dark energy models [TB] . 

To explore Noether symmetry, it is first required to express the action in canonical form, which is possible by 
introducing auxiliary variable in higher order theory of gravity. This has been attempted by several authors |17] 
for F{R) theory of gravity, treating R — 6(^ + f^- + ^) = as a constraint in the Robertson- Walker minisu- 
perspace model. In the process, the Lagrangian is spanned by a set of configuration space variables (a, d, i?, i?) , 
i.e., R is treated as an auxiliary variable and the result is F(R) oc R^ in the vacuum and matter dominated 
era. Such a form of F{R) produces a oc instead of the standard \/i in the radiation dominated era [T5]. This 
creates problem in explaining BBN, since Universe expands at a much faster rate than the standard Friedmann 
model. Further, it has been shown that [18] under appropriate choice of variable hij = = z , required for 
canonical quantization of higher order theory of gravity in the Robertson- Walker metric |19] . z becomes cyclic 
for F{R) =R^ . It means R^ has inbuilt Noether symmetry for F{R) theory in vacuum or in matter dominated 
era. Now, scalar field appears naturally along with higher order curvature invariant terms in the weak energy 
limit of string theories (2^, 4 -dimensional Brane world effective action [2l], supergravity theories [22] and in the 
Chern-Simons gravitational theories [33]. Therefore, adding more degrees of freedom in the form of a scalar field, 
z is no longer cyclic for F{R) = R^ and situation might improve, yielding a better form of F{R) to explain 
presently available cosmological data. Additionally, the symmetry already obtained for F{R) = f{4))R + , in 
the presence of a scalar field |15| must be reproduced, if canonization treating R as an auxiliary variable is correct. 

Higher order theory of gravity may be expressed in canonical form choosing other auxiliary variables too, different 
from R. All lead to correct classical field equations but different quantum dynamics. Horowitz [23] suggested 
to choose auxiliary variable as the derivative of the action with respect to the highest derivative present in the 
action, i.e., as Q = ■ This suggestion is flawed, since in that case auxiliary variable may be introduced even in 
linear theory as shown by Pollock [25] leading to completely wrong quantum dynamics. Later, in a series of works 
[TO] it was suggested that auxiliary variable should be introduced following the prescription of Horowitz [53] only 
after removing appropriate total derivative terms from the action. Following this technique, we found Noether 
symmetry for F{R) = f{4>)R + [3R^ , R being non-minimally coupled to a scalar field [TF. Thus our intention is 
to see if same result [15] is reproduced taking R as an auxiliary variable, and if Noether symmetry at all exists 
for other form of action. 

Thus in the following section 2, we take up a theory of gravity consisting of scalar field being nonminimally coupled 
with F{R) , follow the usual technique of finding Noether symmetry [17] to end up with the result that Noether 
symmetry for F^nn ^ is obscure. In section 3, following the same technique |17| . for a nonminimally coupled 
scalar-tensor theory of gravity along with F(R) term we obtain the same result that Fjm 7^ is impossible. 
This is surprising, since, as already mentioned, we have explored Noether symmetry for the same action, with 
F{R) = f{(j))R+(3R^ and the Brans-Dicke coupling parameter uj{4)) = 1 , in the conformally fiat Robertson- Walker 
metric [15] . In doing so, we have used an auxiliary variable Q , different from i?, and the Lagrangian was spanned 
by the set of configuration space variables (a, Q, </), d, Q, <^), instead. Thus, the absence of Noether symmetry in 
both the cases automatically raises doubt about the technique of treating -R— 6(^-f^-|-^) = 0, as a constraint 
of the theory. To check if indeed the technique is flawed, in section 4, we take up a set of anisotropic models, which 
also increases the configuration space variables (a, 6, i?, ct, 6, i?), to explore Noether symmetry of F{R) theory in 
vacuum. Here again we find that F[R) does not admit any nonlinear form. Thus we conclude that symmetry 
in isotropic space-time obtained for F(R) = R^ in vacuum and matter dominated era, was an accident, since it 
makes z = cyclic, but the technique [17] is flawed, and it is not possible to explore Noether symmetry of F(R) 
theory of gravity, in general. 

The next question that automatically arises, "is it possible to find an integral of motion in general for F(R) theory 
of gravity being coupled with a scalar field, including a linear term in addition?" . The answer is yes and already 
appears in the literature [26] . where a general conserved current in connection with a general higher order theory 
of gravity coupled with a dilatonic scalar has been explored in a metric independent way. In section 5, we briefly 
discuss the technique of finding the same and as an example, show how to use such conserved current to find 



solutions in the context of cosmology. We have given a simple solution of the scale factor directly for the action 
containing F{R) — being added to a linear term (/(0)_R) along with a scalar field 0, in the Robertson- Walker 
metric, both in radiation and matter dominated era. It produces Ftiedmann type solution in the radiation era and 
a coasting solution in the matter era, which fits Snia data perfectly although early deceleration remains absent 
creating problem in structure formation. Finally, we end up with conclusions. 



2 F{R), non-minimally coupled with a scalar field. 



Here, we start with the following action, in which F{R) is non-minimally coupled to a scalar field as. 



S 



d X\J —g 



In the Robertson- Walker line element 



1 — kr^ 



r^dQ^ -t- sin^ 



(1) 



(2) 



„ , a a k 
R = 6{ - + — + — 



As already mentioned, the technique [17^ under discussion is to treat R^^(^^ + ^ + -^ 
to introduce it in the action via a Lagrange multiplier A as. 



(3) 



= as a constraint and 



S 



d X\J —g 



h{(l))F{R) - X\R~6 



a a k 



and to vary the action with respect to i?, which yields, 



(4) 



A = h{<i>)F^R{R). 



(5) 



Next step is to substitute the form of A into the above action. After integration by parts to get rid of the 
appropriate surface term, the action is automatically expressed in the following non-degenerate canonical form, 
that can be checked by calculating the determinant of the Hessian % — g^.J^. | 7^ 0. Thus, 



S 



ha^iF - RFr) - ehad^Fn - 6ha^dRF,,R - 6h'a^d^F, ^ 6khaF, ^ (|^^ - ^ 



dt. 



(6) 



We now move straight into the process of finding a non-linear form of F{R), along with w{(j)) and V{(j)), by 
imposing Noether symmetry, £xL — XL — 0, where, X is the vector field and L is the point Lagrangian. So, as 
usual, let us equate the coefficients of d^, R^, 0^, d(j), (j)R, Rd and the terms independent of them, separately to 
zero to get the following seven equations respectively. 



ah{P + af3,a)F^RR + {ah + 2aa,ah + ajh' + a^-f,ah')F,R = 0. (7) 

6ha^FMRa,R = 0. (8) 

3{wa - ih'F^Ra') + a{w'j + 2wj') = 0. (9) 

a[h'P + hl3']F^RR + [h'{2a + aa,a) + 2ha' + a{h"j + h'i)]F^R - \a^w^,a = 0. (10) 

eha'F^RR + Gh'a^RF^R = waj^R. (11) 

PhaF^RRR + [2ha + haa^a + ha(i^R + h'aj]F^RR + {2ha.R + h'a'y^R)F^R = 0. (12) 

/3/ia[6fc - a'^R]F^RR + a^iF ~ RFM){iah + jah') + ekF^iah + -/ah') - a^{3aV + -/aV) = 0. (13) 



In the above set of equations and everywhere else dash ( ' ) represents derivative with respect to (p . In the following 
subsection, we try to explore Noether symmetry from the above set of equations (7) through (13). 



2.1 Solutions. 



Since, we are only interested in a non-linear form of F{R), i.e., Fj^ 7^ 0, so equation (8) implies a^n = 0, i.e., 
a = a{a, </>) . Hence, equation (11) leads to two distinct cases, which we explore underneath in the subsubsections 
(2.1.1) and (2.1.2). 



2.1.1 Case 1. 

The first case that arises in view of equation (11) is, 

^ aw 

which can be integrated to yield 

T-, o-w . Qha' , ^ , 



aw 



(14) 



where, m is clearly a constant, since F^n , is a function of R only. Hence, substituting this form of 7 and 7' in 
equation (9) we obtain, 

2hF^R{w'a' - 2wa") = w'^a - 2m{2wha" + 2wh' a - w'ha'), {i) 
which holds only under the following conditions, 

2wa" = w'a' => a'^ = Qi(a)^w, (ii) 

w'^a = 2m{2wha" + 2wh'a' - w'ha'). {in) 

These pair of equations relate w , a and h . In fact using (m) , equation (iii) takes the form (Note, a' = => w = 0, 
for which (p turns out to be a Lagrange multiplier, and so we omit this case), 

wa = Amh'a' . {iv) 

Finally, equations (m) and {iv) imply, 

aa' 



4mh'. 



iv) 



Since, the right hand side of {v) is a function of cj) only, so a must have the separable form, a = ai{a)a2{(l>) , and 
thus, 



0:2 = 8mh and w — o!^ . 
Now using equations (14) and (15), one can express equation (7) as. 



(15) 



(C + <a)p3FRR + ■^ai,aa2F,R^ 



a 2 



a2FR = 0. (16) 
where, (3 = ^{a, (/))/33(i?) . Equation (16) leads to three possibilities, which we study in the following three sub-cases. 
Subccise-1 

PaF^RR = nF^R, a ^ a{a), na«),a + Q!iQ!2 = 0. 

Substituting the above condition along with the forms of 7 and 7' from equation (14) in Noether equation (10) 
we arrive at, 



6hh'ai ^ Gh'^ai Qhh'aia'2 



+ 



/Qhh'ai 6/i'^ai 6hh'aia'2 



na{h^)' + 2ai{a2h)' + aia'2h - ml —j—^ + , . ... 

\ a2 (as) 



F^R — aia'2hm = 



A 



So, cither = 0, or, m = 0, yielding a2 = 0. Both the cases imply w = in view of equation (15), turning out 
to be a Lagrange multiplier, which is not our concern. Thus we switch over to the following subcase. 



Subccise-2 



C + a^,a =0, ai = 0, 
which is trivial, or, 

C + cf^o =0, a2 = and aai^a + 2ai = 0, 
which is again trivial. So the following subcase is only of importance. 
Subccise-3, 

q being a constant, along with the coefficients, 



(17) 



g(C + ata) + -^cti,aCt2 = and aai„ + 2ai = 0. 
2m 



(18) 



I.e., 



ai = a 



mqqi 



^2 and a^{l3i + aPi^a) = q2 ^ l3i = - 



Q2 

2a^' 



(19) 



In the above we have further separated ^{a, cj)) as ^(a, (j)) = /3i(a)/32(0) . If we now take up equation (10), it leads 
to, 



Qa'hh" , f ha'\' 

qa{h^y + +6h' {— > 

w \ w 



+ 



„ , Qmhh a „ , , f ha 

(2a + aa o + 2a')h Qmh' 

w \ w 



a'h 



Fr + ma h 



a 



= 0. 



(20) 



Since coefficients must vanish separately, so, either a' = =^ w = 0, which is trivial, as already discussed, or, 
a ^ a{a) . But then, to satisfy equation (18), ai = 0, i.e., a = 0, which again leads to a trivial situation. Thus, 
Noether symmetry for F^rr ^ , remains obscure. 



2.1.2 Case 2. 

Equation (11) leads to yet another situation in which, a' = = 7,_r, i.e., a = a(a) and 7 = 7(0, 0). Thus 
equation (9) yields, 



3-+7— +27' = 0. 

a w 



Now, separating ^ as, ^ = ^{a,(j))P3{R) , equation (7) leads to, 



(21) 



O F^RR 

P3-k — = ni 



a + 2aa,a + ^^(t + "^Ta) 



Fr 

while equation (12) yields. 



a(C + <,a) 



(22) 



F^RRR ^ „ 2(a + aa,a + a^j) 
.—f^ r P3,R = ni = . 



^,RR 



a^ 



(23) 



Left hand side of the above pair of equations are essentially the same, and so we have taken the same separation 
constant ni for the both. Equation (10) now takes the following form, 



F.B^ - a'^l.a = 0. 

6 



{hO'l33F,B.R + (2- + a,a)h' + 2-h + {^h'y 
\_ a a 

At this end, plugging in equation (22) in the above equation, we arrive at, 

mihO' + (2- + aa)h' + 2-h + i^h'Y = 0; 7 = 7(0)- 
a a 

In this case, since 7 is only a function of 0, so equation (21) implies 

w' 

a — a(a) — ma and 7 h 27' = —3m. 

w 

Thus, equation (23) yields, 
h' 



27- 



■ ni^ + 4m = 0, 



which clearly indicates that ^ — £,{(!))■ Further, equation (22) leads to, 
h' 

j-r + ni£ + 3m = 0. 
n 

Last pair of equations (27) and (28) together lead to, 
h' ^ 2m 

= -TO, £ = . 

ni 



7- 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



Thus, ^ turns out to be a pure constant and hence, /? — $/33(-R) turns out to be a function of R only. It can be 
checked that in view of these solutions, condition (25) is trivially satisfied. We are finally left with equation (13), 
which can be expressed as. 



/ 6fc \ 1 

{F--,F,)--i3mV + W')-0. 



(30) 



Apparently, equation (30) is satisfied only if V and h are constants, while k = 0. But then equations (27) and 
(28) are satisfied simultaneously, for m = , which makes a = = /3 and satisfies equation (30) trivially keeping 
the form of F{R) obscure. Summarily, we find that, the Noether equations (7) through (13) can not be solved for 
a nonlinear form of F{R) and so we fail to extract Noether symmetry, at least following the technique of treating 
i?— 6(| + |5- + ^) = to be a constraint of the theory. Situation does not improve under the addition of dark 
and baryonic matter. To be more conclusive, we take up a more general action by adding the Ricci scalar in the 
following section. 



3 F{R) in addition to non-minimally coupled scalar tensor theory of 
gravity. 

Here, we take up even a more general form of action, in which F{R) is added to a non-minimally coupled scalar- 
tensor theory of gravity, as 



S 



-gd X 



f{^)R + BF{R) 



(31) 



Note that this is the same action for which Noether symmetry was found earlier taking F{R) = R^ a-priori and 
following the standard technique of using an auxiliary variable Q different from R, such that the Lagrangian 
was spanned by the configuration space variables (a, Q, 0, d, Q, </>) [5]. So, we expect to get Noether symmetry 
for the above action (31) at least for F{R) = R^ , following the technique of treating R — 6(^ + ^ + -^) = 



as a constraint of the theory. This is what we pose in this section. Substituting the expression of the Lagrange 
multipher and removing appropriate surface term, the point Lagrangian corresponding to the action (31) in tlie 
Robertson- Walker space time (2) reads, 



L = &{-fo? - f'aa^ + kf)a + B{F - RFR)a^ - 6B{FRa^ + aaRF^RR - kaF^R) + a^{^ - V). (32) 
Noether equation are, 

(a + 2aa,a){BF^R + /) + B{P + ap,a)aF^RR + f'a{j + aj,a) = 0. (33) 

&Ba?a,RF,RR = 0. (34) 

2,{uja - A fa') + 0(70;' + 27'a;) = 0. (35) 

Bl3aFRRR + B{2a + aa,a + a^,R)F^RR + f'a-f,R + 2(/ + BFR)aR = 0. (36) 

6BF,RRa' - ojarn. + Gf'aR = 0. (37) 

U[af + (/ + BF^R)a'] + 6[a,a/' + 7/" + Bp'F^RR]a - cua\a + 6aj' f = 0. (38) 

a[6kf + GkBF^R + 3Ba'^{F - RF^r) - Za^V] + l3[B{-a^R + 6ka)F^RR] + j[6kf'a - a^V] = 0. (39) 

FoUowing subsection is devoted to find solutions of the above set of equations (33) through (39) under the condition 
Frrt^O. 

3.1 Solutions. 

Equation (34) implies a ^ o>{R) , i.e., a = a{a,(j)) as before, since we are interested in non-linear form of F{R). 
Hence equation (37) leads to two distinct cases, which we explore in the following subsubsections (3.1.1) and (3.1.2). 

3.1.1 Case-1 

Since, a^R = 0, so equation (37) may be integrated immediately to yield, 
&Ba' wa 

where, as before, n is a pure constant, since F^r is only a function of R. Substituting the above form of 7 in 
equation (35), we get, 

2B{w'a' - 2wq"){F^r -n)= w'^a - Af'wa', {vi). 
The above equation is satisfied if both sides separately vanish, which is possible under the following conditions, 

2wa" = w'a' a'^ = ai(a)^w, {vii) 

w'^a = Af'wa'. {viii) 
Equations (vii) and {viii) imply, 

-T = 4/'. (ix) 

Hence a must admit the separable form, a = ai (0)0:2 (?!>) , and so, 

(^2 = 8/, and w = a'2- (41) 
Now, using equations (vii) and (viii) , in equation (32) following equation is found. 



+ aia)ap3FRR + BFr 



, „ \ Gf'aai aCe'2 
[ai + 2aai,a)Q2 H ' 



+ /(ar+2aar,>2- '^"^'""^-°"^ =0. (42) 



where, f3 ~ ^(a, (/>)/33(i?) . Equation (42) is satisfied under different conditions which we discuss in the following 

three subcases. 



Subcase I. 

^=1, (43) 
IB{^ + aia)a + B [(ai + 2aai,„)a2 + '-^'"t,'""^ ] = 0' (44) 
/(ai+2aai,„)a2- '^"^'r""^ =0 (45) 
where, I is a constant. Now using equation (41), one can express equation (45) as 

Oi2{oii + 2aa\^a) — 12Bnaai^aCt2 = 0. (46) 

Here, we consider 7^ because 0:2 = imply a; = as may be seen from equation (41). Thus, equation (46) 
is satisfied only under the conditions. 

The above two equations imply ai = , which leads to a trivial case. Hence, we take up the other alternative in 
the following subcase. 

Subcase II. Equation (42) is also satisfied provided, 

PsF,RR = I = const, (47) 
IB{^ + aCa)a + /(ai + 2aai,„)a2 - '^^'"^'""'"'"^ = 0, (48) 
(ai + 2aai,„)a2 + '^'°"J""^ = (49) 
Again, using equation (41) we have from equation (49) 

2ai + 7aai,a = 0. (50) 

Now, in view of equation (47) one gets, jS^F^RRR + IS^^rF^rr = 0. Using this together with equations (40) and 
(41), equation (36) takes the form, 

7ai + 2oai,o = 0. (51) 
Equations (50) and (51) are simultaneously satisfied only for ai =0, which again leads to trivial situation. 
Subcase III. 

Here, we consider that the coefficients of equation (42) vanish separately, i.e., 

fc , c \p n f ,0 \ , 6/'aai,aa2 n ^z- , o ^ QBnf'aai^aa'2 „ 
(C + aCa)P3 = 0, (ai + 2aai,a)Q!2 H ' = 0, f{ai + 2aai,a)a2 ' = 0. 

UJ UJ 

Note that, the last equation is the same as equation (45) of subcase I, which only leads to trivial result. Thus the 
case under consideration does not yield any symmetry for a nonlinear form of F{R) . 

3.1.2 Case 2. 

Equation (37) is also satisfied under the condition, a' = = 7,^, i.e., a = a{a) and 7 = 7(0,0). Therefore, 
equation (33) may be expressed as, 

+72— +27^ = 0, 

071 UJ 



where, we have separated 7(a, (j)) = 7i(a)72((/)) . The above equation is true only if, 

= ci ana 72 h 27, = — ci (52) 

071 OJ 

Ci being a constant. Now, separating /? as, /3 = ^(a, (/>)/33(i?) , equation (31) leads to 

+ aia)a(33FRR + BF,R{a + 2aa,,) + f{a + 2aa,a) + fa{j + aj,a) = 0. (53) 

Note that, in equation (53) the condition P^F^rr = constant, automatically implies Fr = constant, which is not 
our concern, so equation (53) is satisfied under the following conditions, viz.. 

Subcase I. 

PzF^RR = cF,R, (54) 

c(? + <Q)a+ (a + Saa^a) = 0, (55) 

/(a + 2aaa) + /'a(7 + a7,a) = 0, (56) 

c being a constant. Now, using equation (54), equation (39) can be rewritten as 

-Ba^{ca£, + Sa)RFM + 6kB{ca^ + a)F^R + Wa^aF + a{<6kf - Za^V) + 7i72(6fc/'a - a^V) = (57) 

Here, none of a , ^ , 71 and 72 are function of R and so, the above equation is satisfied provided the coefficients 
of different derivatives of F vanish separately. This yields a trivial result, viz., a = /3 = 7 = 0. Hence we turn 
our attention to the following subcase. 

Subcase II. 

Equation (53) is also satisfied under the following set of conditions, viz., 

PzF^RR = c, (58) 
cB{i + ai,a)a + f{a + 2aa,„) + /'a(7 + a7.„) = 0, (59) 
a + 2aa.a = 0. (60) 

But then, using equation (58) along with the condition ^_r = 0, equation (36) is expressed as, 

2a + aa,a = (61) 

equations (60) and (61) together imply that a = 0, which leads to trivial situation as discussed earlier. Thus 
we fail to extract a non-linear form of F(R) , imposing Noether symmetry in the most general form of an action. 
Here again the situation remains unaltered under the addition of dark and baryonic matter. It appears, that the 
technique of using the form of the Ricci scalar R — 6{- + ^ + = as a constraint in the action is wrong, since 
as repeatedly mentioned, noether symmetry for such a theory with F{R) — R^ already exists in the literature 

m- 

4 In search of Noether symmetry of F{R) in anisotropic models. 

In this section we proceed to apply the same technique in some of the anisotropic models, taking into account 
pure F{R) theory of gravity. As in the case of Robertson- Walker metric, we can expect some form of F{R) here 
too. The form may be different, since the constraint is different. For, this purpose, we take up spatially symmetric 
Kantowski-Sachs (K-S), Bianchi-I (B-I) and Bianchi-III (B-III) metric, which can be expressed altogether as. 



ds^ = ~dt^ + a^dr^ + b^[de^ + fld(j)\ (62) 
where, fk = sin9 ^ k = +1 (K-S), fk ^ 9 ^ k = (B - I), fk = sinhO ^ k = -1 (B - III), and for which 



Under variation of the action 



/ 



-gdt, 



(64) 



with respect to R yields, the same old form, viz., A = BF^r. Thus the action after removing total derivative 
terms yields. 



A = 



ah\F - RFr) - 2F^ii{2bab + ab^ - ka) - 2F^RR{h^a + 2ahh)R 



dt + T,, 



(65) 



where, S = 2F^ji{db + 2ab)b, which is at par with the surface term obtained from the standard variational 
technique. Now, in order to find a form of F{R), if we impose Noether symmetry {£xL = XL = 0) , we get 
following equation, viz., 



a[b^{F - RFr) - 2f fl(62 - /c) - AF^RRbbR] + l3[2ab{F - RF^r) - AF^Rdb - AF,RR{bd + ab)R] + 

-F^RR\Rab^ + 2{2bdb + ab"- - ka)\ - 2R{ab')F^RRR - + ^6 + ^i? j [4f ^66 + 2F^RRb''R\- 

^Rj [AF^Ribd + ah) + AFRRabR)] - (^^d +^b+ ^R^ [2F,RR{b''d + 2abb)] = 0. (66) 

So, as usual, equating coefficients of d^ ,b^ ,db, R? ,bR, Rd and others to zero, following seven Noether equations 
are found, viz., 

F.RR 



_ .Mia 

Frr _ Q+26a,i,+2a/3,, 



a(7+257.,,) ■ 
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Fh ~ ^ 6(27+2a7,„+f)7,6) 

ba^R + 2a^,R = 0. 

26a+6^a,5+2a/3+2a6;8,6+2a57,B" 



F.RRR I F,RR 

F,R Fh 



2afe7 

2,3+6a,<i+2a/3,<i+b7,H 



+ 



ba,a+a/: 
067 



0. 



F,RR + 



b'^Ra-2ka+2abR0 
a(6^iJ-2fe)7 



67 
F,R 



■2§f =0. 



b(ba+2al3) 
a{b'-'R-2k)j 



F = 0. 



(67) 

(68) 
(69) 
(70) 
(71) 

(72) 

(73) 



4.1 Solution of Noether equations. 

In order to find a non-linear form of F{R) , we are required to solve the above set of Noether equations containing 
four variables a,/?, 7 and F , under the assumption F^rr ^ 0. For this purpose, let us separate the variables as. 



a = A{a,b)Di{R), p = B{a,b)D2{R), j = C{a,b)D3{R). 



(74) 



In view of the above choices made in (74), equations (67), (68), (69), (70) and (73) take the following form 
respectively. 
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D3F,RR 
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a{C+2bC,b) 

bA,„Di + (B+aB,a+bB^i)D2 
b{2C+2aC,a+bC_b) 



D^F RR + 



bADi^R + 2aBD2,R = 0. 

b{bADi+2aBD2)(RF,R-F)-2kADiF^ii 
aC{b-'R-2k) 



0, 



(75) 

(76) 

(77) 
(78) 
(79) 



1 r\ 



where, rii is a separation constant. In view of (76) and (77), one immediately gets 



D2 = miDi, 
In view of which equation (78) now reads. 



(80) 



bA + 2miaB = 0, 



(81) 



where, mi is another separation constant. Replacing F^jiji by F^n and using equations (80) and (81) in equation 
(79) one gets. 



2kA 



mini 



Fr = 0. 



a{b^R-2k)C 

Since 7^ , so the above equation clearly yields a solution of the curvature scalar as. 



(82) 



2k 



A 



miUiaC 



n2, 



(83) 



where, 712 is yet another separation constant. Thus R turns out to be a constant, which is zero for k = 0. This 
result directly contradicts higher order theory of gravity. The field equation for higher order theory of gravity is. 



F,B^ R^^ - -Fg^u + {F,R )'".^ g^i, - {F,r ).^.^ = 0. 
Trace of which yields. 



(84) 



i?F,H+3(F,;j)'".„-2F = 0. 



(85) 



Thus, R can not be zero or constant, which implies Birkhoff 's theorem that Schwarzschild's solution is the unique 
spherically symmetric vacuum solution does not hold for F{R) theory of gravity. In fact, substituting the above 
result that i? is a constant in equation (79), one obtains 



F,B.R = 0, (86) 

which contradicts the starting assumption that F_r_r 7^ 0. Thus, we observe that Noether symmetry of F{R) 
theory of gravity fails to yield a non-linear form of F{R) in the anisotropic models under consideration. Thus, 
despite rigorous calculation, Noether symmetry of F{R) theory of gravity remains obscure. This is not a problem 
whatsoever, since nature may not allow such additional symmetry. Problem is, Noether symmetry for F(R) — R^ 
already exists [15] and the present technique [TT] fails to reproduce it. Hence, fixing F{R) a-priori and expressing 
the canonical form of the Lagrangian using auxiliary variable Q [12] , Noether symmetry may be realized. On the 
other hand, since there is no way to find a canonical point Lagrangian for F{R) theory of gravity, other than 
the technique under consideration [T^, hence at this end we conclude that Noether symmetry of F{R) theory of 
gravity is obscure and so there is no way to find a form of F{R) , other than choosing it arbitrarily. 



5 A general conserved current for F{R) theory of gravity. 

Noether symmetry not only fixes the form of a parameter but also yields an integral of motion leading to a cyclic 
variable and thus helps in finding some exact solutions of apparently non-tractable field equations. Thus absence 
of Noether symmetry for F{R) gravity makes it practically impossible to solve the field equations even if a form 
of F{R) is put in by hand. In this section we show that the situation is not as bad as it appears, because a 
general (Non-Noether) conserved current exists for higher order theory of gravity [53]. For the purpose, here we 



take up action (31) in addition to a matter Lagrangian Lm, to briefly illustrate the existence of the said integral 
of motion, as 



A = J R + BF{R) - uj{cj)) 0,^ f/)'^ - V{(j)) - Ki„,] d^x. (87) 

Field equations corresponding to the above action under metric variation method are, 



+B 



--T 



Rf + 2uo(t)\^;^ + ~ y (0) = 0. (89) 

The trace of the field equation (88) is, 

Rf - 3/;^^ - -2V- B[R{F,n ) + 3{F,r )^"^„ - 2F] = -^t;. (90) 

Combining equations (89) and (90) it has been shown [15] that under the following condition, 

B[R{F,n ) + S{F,n T ,^ ^ 2F] = | t; + J {j-^j ^ (91) 
the above action (87) carries a conserved current 

^;'; = [(3/" + 2/w)^<^^^];^-0, (92) 
which in the case of homogeneous cosmology, takes the following simplified form. 



V-g(3/'2 + 2/ to) = c, (93) 

where, c is a non-vanishing constant and g is the determinant of the metric. Thus, a general conserved current 
exists for higher order theory of gravity, without even fixing the forms of the parameters involved, which lies 
beyond the scope of Noether symmetry, since it is true for the theory of gravitation in general. 

3 

5.1 Attempting a solution for F{R) = R2 . 

It might still appear impossible to solve the cosmological field equations using the conserved current (93). However, 
it has been expatiated in a recent letter [27], that it may be useful to solve the field equations taking F{R) = R^ . 
Here, just to cite another example we try the same for F{R) = R2 . Remember that F{R) = R^ yields a 
Noether symmetry in the Robertson- Walker space-time only in the absence of a linear term and the dilatonic 
sector. However, here we are presenting a solution corresponding to an action containing a linear term being 
non-minimally coupled to a scalar field in the presence of F{R) = R2 . In the Robertson- Walker space time, we 
are required to solve the four available equations, for a{t),(j){t), f{(j)),uj{(j)) and V{(j)). Thus one has to make yet 
another assumption. To handle the equations comfortably, we therefore split the condition (91) into two, assuming 
V^fW\ as 

B[RF,R +3{F,R )^"^„ - 2^] = f Ti:, (94) 
V^Xf, (95) 



t r\ 
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Figure 1: Distance modulus {M — m) versus redshift z plot of the present model (blue) satisfying coasting solution 
a oct, shows perfect fit with the ACDM model (red). 



where, A is a constant. So finally, in the case of homogeneous cosmology, we need to solve following set of five 
equations. 



B[RF,R+3DF,R-2F] ^ § T/^, 

V = xp, 

i?/ - 30/ - - = 0, 

Rf + 2wn(/) + - V'{(t)) = 0, 

V-.g(3/'2 + 2/c^) = c. 

Now for F{R) = i?i , 

dfm = -[(i? + eR)F,R,R + r^Frrr] = 

and so equation (96) reduces to. 



3_R__3i?2_ 27 i? 



(96) 
(97) 
(98) 
(99) 
(100) 



(101) 
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Vr 



9 i?2 
2i?l 



2i?^ 



2k , 



(102) 



Taking Baryonic and non-Baryonic dark matter in the energy momentum tensor, one can solve the above equation 
(102), both in the radiation and matter dominated era in the flat case fc = 0. In the radiation dominated era, 
taking p = ^p, the trace of the energy momentum tensor Tj^ vanishes and equation (102), gets solved to yield 

a cx Vt. This is a brilliant feature of higher order theory of gravity, that despite being tightly coupled to the 
scalar field it still produces Friedmann type solution in the radiation era keeping standard big-bang-nucleosynthesis 
(BBN) unchanged. Note that we are required to find cj), /{(/)), V{(p) and uj{(f>), in view of the four equations (97) 
through (100), which although we are not attempting, may be found in principle, since explicit form of R is now 
known. In the matter dominated era, on the other end, p = and the trace of the energy momentum tensor 
is simply = p, with a^p — PmQ, Pmo being the present matter density of the Universe. The only solution 

that we could find for equation (102) is a = aoi, provided the condition B = — ^q^"" holds. However, such a 
coasting solution although fits Snia data perfectly as illustrated in figure-I, is ruled out, since it does not admit 
structure formation. In our earlier work 18 , we found that R^ yields a on in the radiation era, which creates 
problem in explaining BBN and its solution in the matter dominated era, although admits early deceleration and 
late time acceleration, does not fit with the Hubble parameter versus red-shift - standard model data. Therefore, 
we have tried here to incorporate a linear term in the action, being coupled with a dilatonic sector. It shows, 
despite improvement in radiation era, it does not admit experimental data in the matter dominated era. Hence, 
i? 2 term is likely to be ruled out. 



6 Summary. 



Both the metric and Palatini variation methods have vices and virtues and it is a matter of taste to adopt one. 
However at the end one requires a form of F{R) , any nonhnear form of which makes field equations extremely 
difficult to solve. Noether symmetry not only fixes a form out of indefinitely many, but also yields an integral 
of motion which makes field equations somewhat tractable. For such purpose one requires a canonical point La- 
grangian. The only technique known so far to express F{R) in such a form is to introduce R — 6(^ + ^ + ^) = 
as a constraint in Robertson- Walker metric, which gives correct field equations, in view of the auxiliary variable 
R. Following this technique, Noether symmetry yields F{R) — i?2 both in vacuum and matter dominated era. 
Although such a form is reasonably good in the early Universe, suffers from the disease that it does not produce 
standard Friedmann radiation era |18j . In anticipation that a linear term in addition could have cured the disease, 
here we have enlarged the configuration space by adding a scalar field. Rigorous and detailed calculation shows 
that Noether symmetry remains obscure for Fi?_R 7^ 0. We have also attempted to find Noether symmetry for 
pure F{R) gravity in the background of some anisotropic models and ended up with the same result. Nature 
might not allow additional symmetry, but the problem is that Noether symmetry for a scalar-tensor theory of 
gravity being coupled with F{R) = R^ , already exists in the literature [15], while the present technique fails to 
reproduce the same. This raises doubt about the technique |17) under consideration. 

Obviously, one should ask "how come the technique then yields F{R) = R^ in vacuum and matter dominated 
era?" . The answer has been given in a recent work |18[ , where it has been shown that under appropriate choice 
of co-ordinate hij — = z , an action containing R^ in vacuum and in the matter dominated era, makes the 
variable z cyclic. Thus, the symmetry obtained in Robertson- Walker space-time yielding F{R) = R^ , is due to 
the very speciality of the isotropic and homogeneous space-time and so is a mere accident. Since there is no way 
to find a canonical point Lagrangian for F{R) theory of gravity, other than the one under consideration, so we 
conclude that such a symmetry is obscure for a general F(R) gravity. 

It therefore appears that there is practically no way to handle F{R) theory of gravity. However, the situation 
is not as bad, since a scalar-tensor theory of gravity in the presence of F{R) , indeed carries a general conserved 
current [26l [27], which is independent of the metric and also on the forms of parameters involved, viz., F{R), 
f{4>), ui{(j)) and V{4>). Here in section (5), we have briefiy presented the conserved current and cited an example 
taking F(R) = R^ , to understand its use in finding solutions for higher order curvature invariant terms. The 
presence of a linear term in the action modifies the solution in the radiation era to the standard Friedmann type 
(a cx y/t). Thus BBN remains unaltered, which is of-course a brilliant result. On the other hand matter dominated 
era admits only a coasting solution in the form a act, which although fits SNIa data perfectly, does not allow a 
transition from early deceleration and so i?2 is not a good choice to explain late time cosmological observations. 
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